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We consider the adiabatic pumping of charge through a mesoscopic one dimensional wire in the 
QQ , presence of electron-electron interactions. A model of static potential in-between two-delta potentials 

. is used to obtain exactly the scaterring matrix elements, which are renormalized by the interactions. 

' Two periodic drives, shifted one from another, are applied at two locations of the wire in order 

^SJ i to drive a current through it in the absence of bias. Analytical expressions are obtained for the 

pumped charge, current noise, and Fano factor in different regimes. This allows to explore pumping 
for the whole parameter range of pumping strengths. We show that, working close to a resonance 
is necessary to have a comfortable window of pumping amplitudes where charge quantization is 
close to the optimum value: a single electron charge is transferred in one cycle. Interactions can 
improve the situation, the charge Q is closer to one electron charge and noise is reduced, following 
a Q (e — Q) behavior, reminiscent of the reduction of noise in quantum wires by T (1 — T), where T 
is the electron transmission coefficient. For large pumping amplitudes, this charge vanishes, noise 
also decreases but slower than the charge. 

r-| , PACS numbers: 73.23.-b, 72.70.-fm, 71.10.Pm, 05.60.Gg 
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I. INTRODUCTION 



. The suggestion that electrons can be supplied one by one by a mesoscopic circuit has been proposed over two 
decades agc^. Instead of applying a constant bias voltage to the system, it is possible to supply a.c. gate voltages 
^ , which perturb the system periodically. Under certain conditions, the charge transferred from one lead to the other, 
''O ' during one period, can be almost quantized. Adiabatic pumping of electrons could in principle be used in future 
^ ] nanoelectronics schemes based on single electron transfer, and it also has applications to quantum information physics. 
Over the years, theoretical approaches to this adiabatic pumping based on scattering theory have become available^ii^ii^. 
These situations typically describe mesoscopic systems which are large enough, or sufficiently well connected to leads 
that electronic interactions (charging effects for instance) can be discarded. Scattering theory has been applied^ to 
calculate both the charge and noise in systems in the absence of electron-electron interactions. 

On the experimental scene. Coulomb blockade effects have been successfully exploited to achieve pumping with 
isolated quantum dots^. To our knowledge, pumping experiments which are not entirely based on Coulomb blockade, 
^C^' [ where the shape of the electron wave functions is modified in an adiabatic drive are rather scarce. A recent study^ 
fS| ■ has dealt with the transport through an open quantum dot where such interactions are minimized. 

Besides Coulomb blockade physics, the effect of electron-electron interactions in conductors with reduced dimension- 
ality have been discussed by several authors. The case of strong interactions in a one dimensional quantum wire was 
presented in Ref. 0, using Luttinger liquid theory. Alternatively, Ref. |8i discussed the opposite limit, where the effect 
, of weak interactions can be included in a scattering formulation of pumping using renormalized transmission/reflection 
^ ■ amplitudes^. However, the results for the pumped charge remain mostly numerical in this work. 
. ^ The conditions under which pumping amplitude and interactions must be tuned to achieve quantized pumping are 

not obvious. Many physical parameters enter this problem, such as the amplitude of the pumping potentials, the 
H phase difference between these, the possibility of a constant offset on these potentials, the overall conductance of the 
. . 1 unperturbed structure, and to what extent the strength of electron-electron interactions play a role. Analytical results 
on this issues are highly desirable, as well as information about the noise. 

With regard to the experiment of Ref. 0, there is clearly a need for further understanding the role of weak interactions 
in such mesoscopic systems in the presence of pumping. The purpose of the present work goes in this direction, in 
the sense that we provide analytical expressions for the pumped charge and the noise for a one dimensional wire 
in the presence of interactions. This allows us to explore all pumping regimes^^ (weak to strong pumping) and to 
determine in which manner and to what extent the pumped charge can help to achieve single electron transfer. Besides 
addressing the question of the ideal conditions for good charge quantization, we shall establish relationships between 
charge and noise in different regimes. For concreteness, a two-delta potentials model will be used and interactions 
will be added on top of it. 
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II. PUMPED CHARGE AND NOISE 
A. Adiabatic pumping in non interacting systems 

Here, we recall the formula which was established for the charge transferred during the single period of an adiabatic 
pumping cycle through a quasi onc-dimensional system. The system is in general described by a potential V{x) con- 
taining two internal parameters which are modulated periodically. The time dependence is assumed to be sufficiently 
slow so that, although the scattering matrix depends on time, its variations are minute when an electron is scattered 
in the mesoscopic wire. 

At finite temperature, the pumped current readsii 



dt / f{E)Tr{S\E,t)<j,S{E,t)-l} 



dE 
2^' 



(1) 



where S{E,t) is the Wigner transform of the scattering matrix S{t,t'), and to, a good approximation the scattering 
matrix for the problem with the potential frozen, az is the ususal Pauli matrix, / the identity matrix and / the 
Fermi-Dirac function. 



/oc 
e''^^'-''^S{t,t')dt'. 
-OO 



(2) 



The pumping potential will generate sidebands at E±huj and we assume that temperature is much smaller than uj, 
i.e. ksT <^ huj, so that we can approximate the Fermi function by a step function. In fact temperature dependence 
occurs in two places in this problem. First in the Fermi function and second, the scattering matrices elements depend 
on the temperature because of the renormalization due to the interactions (see next chapter) . Formulas for averaged 
current and zero frequency noise can be carried out using results of the literature in the "zero temperature" formalism, 
except that the scattering matrix elements are in fact temperature dependent. 

The pumped charge reduces to a time integral over a pumping cycle^ii^ 
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dX 



dt 



where Su {i — 1, 2) are the elements of the scattering matrix s{E): 



dt. 



(3) 



s{E) = e"^ 



(4) 



where 4> is the phase accumulated in a transmission event and 6 is the phase characterizing the asymetry between the 
refiection from the left-hand-side and from the right-hand-side of the potential. Conservation of probabilities imposes 
R + T = 1. We assume the quantities Vt, and (j) to be functions of the Fermi energy Ep and of the external 
time- varying parameters X{t) and Y{t). 



B. Inclusion of weak interactions 

In the case of weak interactions, the transmission and the reflection amplitudes Si2 and sn can be calculated 
in the presence of Coulomb interaction via a renormalization procedure^. High energy scales above a given cutoff 
are eliminated. The high energy cutoff is lowered progressively. The renormalization has to be stopped when the 



temperature becomes comparable to this cutoff. Finally, if Sj^^'' and s^i denote respectively, the transmission and 
reflexion coefficient without interactions, si2 and sn can be expressed in the formSi: 



Sl2 
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where I = ksQ/W] 8 is the temperature, ks the Boltzmann constant and W the original bandwidth, a is a negative 
exponent related to the strength of the screened Coulomb interaction potential^. Specifically: 

Vc{2kF) - K(0) 

a = 7, . (') 

ZTTVp 

with Vc{q) the Fourier transform of the screened Coulomb potential at g = 2kp and g = 0, respectively, T4(0) is 
finite due to screening, a — Q corresponds to the absence of electron-electron interactions. Tq — js^^"*!^ represents 
the conductance of the wire in units of /h. ^From now on, Q will denote the pumped charge with interactions and 
Qo without interactions. The integrand of Eq. ([3]) is therefore modified by the presence of the interactions. Thus, 
temperature dependence occurs through the renormalization of the S'-matrix. We shall be interested in the regime 
where temperature is much lower than pumping frequency, fc^B ^ tj, as said before, but renormalization of 5-matrix 
should not be too severe so that the renormalization of S'-matrix still makes sense. At very low temperature, all 
barriers become almost opaque and bosonization is required^, so typically > 10^^, that is ks e > VFlQi/". For 
example, for nanotubes having a around —0.3, this gives Itlo 3> fc^G ^ 10~'^VF. 



C. Two-delta potentials model and pumped charge 

Consider the pumping charge Q transferred during a single period through a ID chain of an arbitrary potential 
shape. Let two parameters of the system be modulated periodically. The single particle Hamiltonian reads 

H = ^ + V{x) + Vp{x,t), (8) 

where Vp{x,t) is the time-dependent perturbation part of the arbitrary potential and has (5— like potentials form 

Vp{x,t) = 2kFX{t)5{x ~ X,) + 2kFY{t)5{x - Xf), (9) 

with the amplitudes X{t) = Vi{t)/2kF , Y{t) = Vi{t)/2kF, where Vi{t) and Vf{t) have periodic time evolution with 
the same period to = 27r/a; and kp denotes the Fermi wave-vector,. V{x) is the static potential in between two 6 
potentials. Below, superscript indexes (0) indicate non-interacting systems. 
Using the known rclationi^, 



s 



(0) 

af3 



^Sa0 + 2ikFG°{xa,xp), (10) 



where G'^lxonX^) are the usual real space retarded Green's functions and the fact that the functional derivative of 
the Green's function /dV{x) can be written as the product of two Green's functions^ 

^^lv{x,T^ -G"(x.,x,-)G°(.,,x,). (11) 
For the first bracket of Eq. ([3]), we get 

dx + dx 12 - 2ikF ' ^ ' 

where we used the condition s^i^'s^'^'' * -I- s^i2 * — 1> ■^12^ a-nd s^fi are the bare transmission and reflection amplitudes 
from the disordered system, without taking account the electron-electron interactions. For the second bracket of Eq. 
we get 

,(0) , , aJo) , , .(0) 



5sii (o)» 5si2 (0), _ / (0) (0)* (0)*.. (o)A _ |To, 

^^11 +^^12 -^1^12*11 +^12 (1 + "22)J-^- (13) 



Here, we used the current conservation requirement 521^^2'^^ * -t- 3^12 *S22^ = 0. s'^2 the reflection amplitude from the 
right of the scaterrer and can be presented as 

4"2^=2...^-l. (14) 
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Substituting Eq. ([T2|l and (fT3|) into Eq. Q, using the identity p4)) . we finally arrive at 



Qo = - 



27r 



(15) 



Similar expression for Qo, can be found in Ref. Il4 

Using the same method for Q, we obtain, in the presence of interactions 
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For numerical simulations, we specialize to the case where the static potential V{x) = in Eq. ([8]) and the time- 
dependent perturbations are two 5-like potentials separated by a distance 2a. The expressions of 5^2'' ^^.d s^^^^ are 
needed, when the single particle Hamiltonian reads 



H ^ -h^k^ /2m + Vp{x,t). 
The elements of the 5- matrix in the absence of electron-electron interactions, are given by 

(F - X)sm(2fcFa) - i{2'XYsin{2kFa) + (X + Y)cos{2kFa)] 



SO) 
'12 



D 



1 



with 



D = (^l-2XYsin^{2kFa) + i[X + Y + XYsin{4:kFa)]^, 



and X — X/2, same for Y. s'^^ — s^2^ and 532^ is obtained by replacing F by X in s^^^ . 



D. Noise 



(17) 

(18) 
(19) 

(20) 



Ref. [ij studied adiabatic quantum pumping in the context of scattering theory. Their goal was to derive under 
what conditions pumping could be achieved optimally, in a noiseless manner, with the assumption that the pumping 
frequency is small compared to the temperature. This enabled the authors to derive expressions not only for the 
pumped charge per cycle, but also for the pumped noise, the current-current correlation function, averaged over a 
time which is long compared to the period of the adiabatic drive, at zero frequency. Specifically, the noise is defined 
from the current-current time correlator: 



S{t,t') - ^{SI{t)SI{t') + SI{t')Sm), 



(21) 



with 61 — /—(/). This correlator is then averaged over ng periods of the pumping drive with no large, and it is 
taken at zero frequency by performing an integral over the remaining time argument. Setting tq — tt-q 27r/a;, 



s(si = 0) 



27rrio 



dt 



dt'S{t,t') 



(22) 



Ref. "5 extended the results of Ref. 14 to the case where this limiting assumption is relaxed, yielding a complete 
description of the quantum statistical properties of an adiabatic quantum pump, albeit restricted to small pumping 
amplitudes. Results made use of the generalized emissivity matrix. These results were generalized to arbitrary 
pumping amplitudes by Ref. [ill. Here, our goal is to address the question whether electron-electron interactions 
affect the pumping noise and how. Following formula (14) of Ref. [ill and applying their Eq. (15) without assuming 
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that the pumping amplitudes X(t) are small, it is possible to put the zero frequency noise S for arbitrary pumping 
amplitudes into the form 

1 „2 pTo poo poo poo 

^(^ = 0) = TTTT^- / dt d{t'-t) /(-ei) / /(es) 



(2^S)2 To 7o 

X Tr 



— OQ J —OO J —OO 



where s(ei,f) is the 2x2 5'-matrix, for an incoming wave at energy ei + ep and value of the pumping parameter 
X{t). tq is a time which is much larger than the period t. t\ — f?k\l2m — ep, where ki is a wave vector and ep is 
the Fermi energy ep = h'^kp/2m. / is the Fermi-Dirac function. Here, s(ei,<) is taken for fixed t^. Since we work at 
temperature much smaller than hw/kB, as explained before, we can set /(— ei) ~ 1 for ei > and otherwise. We 
set e'2 = —£2 so that both ei and e'2 will be positive. M is defined as 

M(ei,e^,t) = s(ei,i)t(T,s(-e^,i). (24) 

S now reads 

S = TTTvi^^^r.^oo - dt' dej de'2Tr\M{ei,e'2,t)M{ei,e'2,ty-l\e-' h . (25) 

{2nhy Jo roJ_oo Jo Jo ^ J 

Now, for large pumping amplitudes, the above formula needs to be rearranged, using the fact that X{t) is a periodic 
function of period 27r/w. Note that the dependence of M on ei and £3 prevents the direct use of fast Fourier 
transform. Nevertheless, we can use the fact that, for given values of £1 and e'2, M{t) and M(t') are periodic functions 
of t. Switching to Fourier transform 

2tv_ 

Mn{eue'2) = ^ Mie,,e'2,t)e-'"^Ut, (26) 







M{e,,e'2,t) - J2 M„iei,e'2)e'^^K (27) 

n— — oo 

Performing the trace, we arrive at 

■ dej ^ (|Mi,i„|2 + |Mi,2„P + |M24u|' + |M2,2„|'-25„,o)<5(^^-nL^), (28) 



s 



71— — 00 



where 6n,o is 1 if n = and zero otherwise and Mijn is the element of matrix of Af„(£i,£2), where energy 

dependences have been omitted to ease the notations. 

When hu! is much smaller than ep, formula ((28)) can be simplified. In this case, M{ei, e'2,t) will be different from 
M(0, 0,t) only when ei or e'2 are a non neghgible fraction of ep. This occurs because M(0,0,t) corresponds to 
matrix M for incident wave and outgoing wave at energy ep. We denote by eip typical energies of the order ep. 
eip will correspond to n of the order [eip/hw), which is very large. The Fourier transform M„(£i,£2) will decrease 
exponentially with n for large n. Thus, we can neglect the dependence on £1 and e'2 and replace them by zero, which 
amounts to replacing the energies by £i^, except in the argument of the 5 function. Under these conditions, we have 



S ~ e'^ 
27r 



Ln>l 



n ( |Mi,i„P + |Mi,2„P + |M2,1„P + |M2,2nP 



(29) 



With our form of the S'-matrix, this formula is equivalent to Eqs. (24a), (24b) and (24c) of Ref. 15;, apart from am 
overall factor 2. See Appendix C for details. For numerical simulations however, we did not make this simplification 
and kept the dependence on ei and e'2 of Eq. 



III. DISCUSSION OF PHYSICAL RESULTS 

We now illustrate these formulas by computing the charge and noise, in the case of two-delta potentials model. The 
two parameters of the drive X and Y (Eq.(9)) are chosen to vary periodically according to: 

X = Xo + 7]cos{ujt), (30) 
Y = Yo + 7jcosiujt-ip), (31) 
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where Xq is a constant offset potential and Lp a phase difference. Note that X, Xq and i] are all dimensionless, see 
Eq. Q. To ensure maximal pumping, we shall specializei?^ to ip = 7r/2. 



A. Zero offset 
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FIG. 1: Q, charge with interactions, (solid line), Qo, charge without interactions, (dashed line), both in units of e. S, noise with 
interactions, (dotted line) and So without, (dashed dotted line), multiplied by 2tv/uj, in units of e'^, vs. rj. Essential parameters 
are Xq = 0, no offset, — 0.3, kpa — 0.5 and hco — lO^^ep. 



0.09 — 

0.08 - 

0.06 - 

0.04 - 




0.2 0.4. „ 0.6. 0.8. 1. 



FIG. 2: Q, charge with interactions, (pluses), and Qo I (upper dashed line), both in units of e vs. for rj = 0.3, Xq = 0, 
kpa = 0.5 and tuj = lO^^ep- 



First, the case without offset Xq = Iq = is studied. To look at the influence of interactions, we plot in Fig. 1 the 
pumped charge in units of e, with interactions and without, versus the amplitude of the drive rj, for an interaction 
parameter = 0.3 (moderate electron-electron interactions). There are three regimes: weak pumping, 77 <g; 1, 
intermediate pumping, 77 of order 1 and large pumping amplitudes, 77 ^ 1. The current noise times 2t:/u;, in the limit 
of small uj, is plotted together on the same figure in units of e^. Analytically, for 77 ^ 1, Q reads 

Q = ^sm(4fcFa)/-2"r;2. (32) 

As noted in Ref. 8|, in the weak pumping regime, charge Q is larger with interactions by a factor see Fig. 

2. Results of Ref. [j for the noise, valid for weak pumping and no interactions can be adapted in a straightforward 
fashion to the case with interactions. We find the following formula for the noise for weak pumping. 

S = eH-^W — - (33) 
' 27r ^ ' 
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FIG. 3: Q, charge with interactions, (pluses), and Qo l^" , (lower dashed line), both in units of e, vs. for 77 — 15, Xg = 0, 
kpa — 0.5 and huj — 10~^eF- In this regime, Q is approximately larger than Qo by a factor 
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FIG. 4: (27r/t^)»725', noise with interactions scaled by ri^ and by the period, (+) and {2-k / u)) rf SqI^"" , lower solid line, for 
moderate rj = 15. Also illustrated by the top two curves, is the very large pumping amplitude regime, {2-k/u)) r]'^ S (x) and 
{2-K / uj) rf So (*), for very large r); rj = 100. The top two curves are close to each other. An attempt to fit rfS by l^°'rf'So 
(dashed line) for t] = 100, clearly fails for this regime of pumping amplitudes, kpa — 0.5 and huj — 10~^eF for all cases. 



The noise is thus increased by the same factor as the current. The Fano factor, defined as the ratio S/e{I), is 
4/sm(4fcp'a) and remains independent of the interactions, as long as we remain in the weak pumping regime. This 
corresponds to the very left part of Fig. 1, for 77 smaller than 0.25, typically. 

At intermediate pumping amplitudes, Q reaches a maximum value Qmax which is again larger than its non- 
interacting analog Qomax- This maximum is of the order of the single electron charge, but less than it. Meanwhile, 
the noise decreases. This is a reminder of the reduction of the noise by a factor T{1 — T), where T is the electrons 
transmission coefficient in quantum wires. This explains why the noise exhibits a first maximum around 77 close to 
1, since Q gets closer to one electron charge, noise will decrease. Then, for moderate amplitudes, 77 around 6, charge 
decreases and passes through the value 0.5e, this corresponds then to the second maximum of the noise. 

For large, but not very large pumping amplitudes, typically 77 — 10, Q remains smaller than Qq but behaves in the 
same way, namely as 77"^, as noted in Ref. 10. As a function of the interaction parameter, Q behaves as QqP'^. See 
Fig. 3. For very large pumping amplitudes, (typically order 100 or more), Q becomes practically equal to Qo, Q — Qo 
behaves as 77"'*, see appendix B for details. 

For the noise, we found numerically that S and its analog without interactions, Sq, both decrease as 77^^, much 
slower than the charge. See below for analytical derivations. As concerns now the interaction dependence of the noise, 
S is always smaller than Sq, but for very large 77, S tends towards 5*0. More precisely, for large but still reasonable 77, 
of the order 10 typically, S is almost equal to Sol^°', whereas for very large 77, of the order 100 or more, S and Sq are 
practically the same. This is not surprising since Q and Qq are then also practically equal in the end. 

This dependence on the interaction parameter is shown in fig. 4. We have to plot {2tt / uj) Srf , vs. Z^", but the 
overall factor 2-k /uj is unimportant; the product Srf can be compared to both P°'SQri^ and to 5o77^ for 77 = 15 and 
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r] = 200. For rj = 15, we see that rj'^S is fairly well approximated by tj'^SqI'^"'. On the contrary, for rj = 200, such a fit 
fails and instead, r]-^S is almost equal to 77^ 5*0- 

The results at large rj can be derived from analytical formula for the charge and noise. An expansion for large rj is 
performed, X and Y behave as rj, except at particular points where X or y are zero. See Appendix A for details. 



B. Non-zero offset 



We now turn to the case where Xq is non-zero, which enables to have regions where Q is almost quantized. There 
are basically three cases, according to the value of kpa. 

The first case corresponds to kpa — nir /2 (rigorously), where n is an integer. In this case, it is impossible to pump 
anything. The reason is given below. The second case corresponds to the case where kpa is small but non-zero, 0.1 
typically. We first describe the behavior, then give numerical illustrations and last provide analytical justifications. 
In this case, the charge is almost zero up to — A"o. It rises quickly around r/ — ^/2Xq and reaches a value close to 
quantized e for a wide range of values of 77. This is the quantized region of "q. The width of this region can be shown 
to scale approximately as {kpa)~^ . After the end of this region, Q and Qo first decrease abruptly and for even larger 
values of 77, decrease slower, as 77"^. The noise in the quantized region and around it seems to be well approximated 
by (5(e — Q), reminiscent of the noise for fermions in narrow quantum wires. However, this does not last when 77 
becomes noticeably out of the region of almost quantized charge, since Q and Qq behave as 77"^, whereas S and Sq 
decay only as 7/~^. 

Fig. 5 shows noise and charge with and without interactions, versus rj for kpa = 0.1. Fig. 6 shows a comparison 
between S and a least square fit of the form C Q{e — Q), in the quantized charge regime, where C is the only adjustable 
parameter. See captions for details. We now turn to analytical justifications of the previous assertions. 

1.4 



1.2 - 




FIG. 5: Q, (solid line), curve with a plateau for 77 between 28 and 43, (the quantized region), Qo, (dashed line). The two curves 
with two spikes located at ri around 28 and at rj around 43 represent {2n/uj)S, (dotted line) and (27r/cj)So, (dashed dotted 
line), vs. 77. Essential parameters are Xq = 20, — 0.3, kpa — 0.1 and fkj = 10~^eF. Electron charge e is set to 1. The solid 
horizontal line of ordinate 1 and the thick vertical line at = X0V2 ~ 28.28 are guides to the eye. 



We now explain why pumping is impossible for kpa = nn/2. Since sin{2kpa) = 0, the scattering matrix now 
depends only on a single parameter, the combination {X + see eqs. (|18II20I) . so we denote by s[°^' its derivative 
with respect to X + Y. Thus, 

^(X + F)dt = 0, (34) 

because the bracket is just (1/2) d/d(X + Y) {\s^i'\^ + IsiaT)- Then, we look at the case when kpa is close to mT/2, 
but different from it. Clearly, when sin(2kpa) is small, for r]\sin{kpa)\ < 1, we will be back to the former case, so 
Qo can start to level noticeably from only for 77 values larger than a critical value 77ci, which is proportionnal to 
l/\sin{2kpa)\, independently of Xq. There is at least another scale, namely Xq. For Xq large, (typically larger than 
10) and ?7 smaller than Xq, the pumping contour is a circle which does not enclose the origin and both transmission 
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FIG. 6; S (solid line), in units of e^, and best fit of the form C Q{e — Q) (dashed line) in the interval 20 < 77 < 45, for the same 
situation as Fig. 5. The fit no longer works in the large pumping regime, for r) between 45 and 80 here. 



^-nd derivatives of the transmission coefficients, ds^i2 /dX are small. Qq will remain very small. Thus, to have 
a significant Qq, one needs rj > max{Xo, l/\sin{2kpa)\) , where max{x,y) is the larger of x and y. For even larger 
Tj, when terms like sin{2kpa)XY dominate over terms linear in X or Y, i.e. (|sm(2fci?a)| 77) ^ 1, it is possible to 
expand in 77"^ and we are back to the large pumping regime where Q decays as ?7~^. So, for \sin{2kpa)\ smaller than 
Xq, there will be a region between max{XQ, Ci/ sin{2kFa)) and C'2/ sin(2kFa), where Ci and C2 are constants, where 
Qq is appreciable. These very qualitative arguments however do not prove that the charge is almost quantized in this 
interval, whose width is of the order \sin{2kF0,)\~^ ■ 

We now turn to the analytical explanation of the Q{e — Q) behaviour in the quantized charge regime. It is useful 
to disentangle the effects of the fluctuations of T and those of the phase 0. We start from (|29p and use the fact that 
for any periodic function /(a;), if /„ denotes the Fourier component at frequency mj, denoting x = ut, 



E 

n>l 



Applying the former equality to Mi_i„, Mi_2n, Af2,in and M2,2n, and denoting by Mii{x) the inverse Fourier 
transform of Mi_i„, we get 

S ^ c^^l I f|Mii(x)-Mii(xOP + |Afi2(a;)-Mi2(x')P + |M2i(a;)-Mi2(a;')l' + |A^22(x)-Af22(x')|' 



27r 



i^--r'-^- (36) 



An 

In terms of 9 and T, using ^ and this reads 



S = 8e 



(^) (/. + /.), (37) 



with 



2tt p2tt 



\T{x) -T{x')\^ dxdx' 

{x - x'Y 47r2 
\g(x) ~ g{x'y^ dx dx' 



h = I I ' \L I2 -TTT^ (38) 

(39) 



{x - x'Y 47r2 

with g{x^ = •\/r(l — T) e*^. We expand \g{x) — g(a;')p and rewrite I2 in the form 

J1 + J2, (40) 

with 

- y^T(x) (1 ~ T{x)) ~ ^T{x') (1 - T{x')) - ' 



Ji 



(x — x') 



dx dx' , , . 



- I I ^nx){i-nx))^nx'){i-Tix') f ' (42) 
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Thus, noise breaks into three parts, one related to the fluctuations of T, another to the fluctuations of \/T{l — T), 
and the last to the fluctuations of 9, in fact, to the variations of the slope d9/dx, since has to vary by 27r in one 
cycle to get appreciable pumped charge. The physical message is that if the fluctuations of T are much smaller than 
the fluctuations of the phase, then noise shows a Q{e — Q) behaviour. Otherwise, fluctuations of T bring an extra 
noise that does not contribute to the pumped charge and overall noise is thus larger than Q(e — Q). We now try to 
establish this more flrmly. 

In this paragraph, we now show that, for Qo < e/2, 

5*0 > 8(|^)Qo(e-Qo). (43) 
The first integral J\ involves solely the fluctuations of T and can be rewritten as 



Ji = ^n|Vni-T)„r (44) 

n=l 

We have a lower bound for Ji by replacing n by 1 in all terms of the sum except the term for n = 0, 

OO — 

Ji > Y.WT{l-T)S -WT{l-T),\\ (45) 

n=0 

Using then the Parseval identity for the first term and using the notation (/) = f{x) dxj^-K for the second, 

Ji > (iVT(l-r)l') - l^^T(\-T))^ = (T) - (r2) - (Vr(l-T))'. (46) 
Now for J2, applying twice the Holder inequality, 

r'^'^ Hr' I I \p-iS{x) _ f,-ie{x')\2 

r"^^ dr' I / , \ Ip-^^(^) - p-«^(^')|2 

>- I ^Vnx')(i-n.'))(v/T(r^)/ 

> (v/T(i-r)) -J^ ^-—^^ . (47) 

For the last double integral, we proceed as before, going again to Fourier transform, isolating the component of order 
and using the Parseval identity, it is larger than (|e*^p) — |(e*^)p. Always, je'"! = 1 and, for reasonable 9{x), we 
can assume symmetry x into —x which implies {sinO) = 0. We can also assume symmetry when x is changed into 
■K — X, which implies {cos6) = 0. We have 



J2 > [^yT{l-T)j . (48) 

Now, for 7i , using the same technique, (going to Fourier transform and isolating the n = component) 

h > (T^) - {Tf. (49) 

Putting everything together 

h + i2> (r2) - (r)2 + (T) - (r2) - ^v'i^(i-r))' + (v^t(i^)' = (r) - (r)^. (50) 

Now, assuming that J^^ ^ if ~ ^' (*^^ circulation of (j) is zero in one cycle), 

e-('-(^£))- 

which implies 
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Moreover, by the Holder inequality 



Now, for any y > 1/2, y{l — y) is a decreasing function of y. Thus, if {T^\ > 1/2 , i.e. Q < e/2, then 



Thus 5 > 8(^)go(e - Qo), which is 

Now we turn to the case of interest, when Qo > e/2, for example in the quantized region. We were not able to 
provide a general analytical proof of Qo(e — Qo) behaviour. We first look at simple limiting cases. Then, we examine 
the particular case of the two delta-potentials model. 

First, in a situtation where T(x) is a constant T, (with T small to have almost charge quantization), then Qo = 
e(l — T) and Sq = 8Co(^)Qo (e — Qo) with Co depending on the shape of 0{x), but always Cq is greater than 1. In 
the case of constant slope dO/dx — 1, Co = 1. Second, in the situation where dd/dx is constant but T{x) arbitrary, 
we have (T {dO/dx)) — (T) and thus again (|43p holds for any Qo, not just for Qo smaller than e/2. 

However, in practise, 9 varies abruptly by tt in the vicinity of resonances. This is different from optimal pumping 
strategies which have been studied before^^ii^. We first give qualitative arguments and then give precise calculations 
for the model studied here. Let us look at the contribution of J2, Eq. to the noise. When x and x' are both close 
to resonances |e~'^(^^ — e"'^^^ ^p/ {x — x'Y behaves as (dO/dx)'^ for (x — x'){d9/dx) < 1. If meanwhile, T{x) does not 
vary too much and assumes the value Ti, the contribution of this region in the plane (x, x') to J2 will be Ti{l — Ti)/A. 
The 1/4 comes from the fact that 9 varies suddenly only by tt and not by 27r at each resonance. Another contribution 
will come from the regions where x is within the resonance but x' just outside or vice versa. Let us take x' outside 
to be precise. Then |e~*^*^^) — e"*''^^ '1=2 and integration on x and x' will give a contribution mainly from x' just 
outside, due to the rapidly decreasing factor (x — x')"^. This will eventually give another factor Ti{l — Ti)/A. Despite 
the non-local character of the integrand in J2 , regions where x and x' are both far from a resonance make very little 
contribution to J2- For Q, one gets Q = e(l — Ti/2). In the simple case where there are only two resonances and 
the TiS are equal, then Q = e(l — Ti) and J2 = Ti(l — Ti). If Ii and Ji, which are related to the fluctuations of T{x) 
are much smaller than J2, then, this leads to So = 8 (a;/27r) Qo(e — Qo)- 

We now test the former very qualitative ideas by analytical calculations on our particular model. A first thing to 
be noted is that the phase (p does not intervene in the noise, which is normal since noise is related to the two-particle 
scattering matrix. On the contrary, it does formally enter the equation for the pumped charge, see Eqs. (Pl). (jl5p and 
(|16p . Nevertheless, the variation of (j) when Lot varies in one period, is always zero so that (p does not play any role. 
This is due to the fact that (p is the phase of D, see Eqs. (fT9|) and ([20]). If the real part of D becomes negative, then, 
its imaginary part cannot be zero and thus, 4> can never be equal to — tt or tt. Moerover, varies by strictly less than 
27r during one period and thus, its circulation is zero, giving no contribution to Q and Qo. 

We now turn to the variations of 9. Charge quantization necessitates 77 larger than \/2Xq and kpa small. We thus 
approximate cos{2kFa) by 1 and set u — sin{2kFa) <C 1. can be written as 

9 — arg{n), (55) 

with 



n = XYu + X + Y + iu{Y - X). (56) 

There exist two values of 77, 771 and 772, which play a particular role. For 771 < 77 < 772, when cut varies by 2tt, Re{n) 
changes twice its sign and 6 varies by 27r. In the model studied here, 771 = Xq\/2 and 772 = Xo\/2 + 2\/2u^^. Outside 
this interval of 77, the increase of 9 when ut varies by 27r is zero, not 27r. The reason is the following. For 7/ < 771, Re{n) 
is always positive and the phase 9 remains confined in an interval contained in [— f , f ]. For 7; > 772, Re{n) changes 
four times its sign but the contour described by n, in the complex plane, as ojt is varies by 27r, does not enclose the 
origin. It can be seen directly, for it is impossible to have Im{n) = and Re{n) < at the same time. Thus, it is not 
surprising that the quantized region corresponds approximately to the interval [771 , , 772] . 

We now examine more precisely the variations of 9 and T. Apart from a small variation around x = tot — — 37r/4, 
9, as a function of x, is essentially flat, except around two points xi and X2. In practise, xi is close to — |- and X2 
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close to 37r/4. Around those two points, 9 increases fastly by almost tt each time. We can thus model the function 6 
by 



— vr, for X < xi — I, 



, ior X ^ XI — I, 

X — Xl \ 

— J 1 j for Xl — Z < a; < a;i + Z, 



TT 

2 

0, for Xi+l<x<X2 — I, 

TT / X — X2\ r , 

— (^1 H J j tor < X2 — I < X < X2 

for X > X2 + I, 



+ 1, 



TT 



(57) 



I being a small distance. As for T, T{x) shows a large peak around x = Xs = — 37r/4, and two smaller peaks, 
practically identical, centered around xi and X2- Away from these values of x, T{x) is almost zero. 

Now, we look at the implications for the charge and noise. When calculating the pumped charge without interaction 
via the formula Jq^ ^ {^~'^ix)) dx, the region around x ~ Xg does not bring much contribution because variations 
of 9 are small here. For the calculation of Qo and Sq, we can ignore the large peak in T and thus model T{x) by 



T{x) 



Ti exp 



(X - Xif 



{x - XiY 
T{x) = otherwise. 



for |x 



< I 



(58) 



T has to be derivable in order to avoid logarithmic divergencies due to the factor l/(x — x')^ in Eqs. ([38]) and ([39 
Then, it is possible to perform analytical calculations which give 



go = e(l-CiTi), 

5*0 = C2[^){e-Qo)[e{l~Ci)+CM]- 



(59) 
(60) 



Detailed calculations, involving the integrals Ii and I2 and the constants Ci, C2 and C3 are given in the appendix. 
C3 is smaller than 1, (approximately 0.58). Even if it is not exactly of the form Qo{e — Qo), when there is good charge 
quantization, i.e. when Qo is not far from 1, So goes as (e — Qo)- The essential thing is that T and 9 vary rapidly 
around certain values of cut. We do not get exactly Qo(e ~ Qo) because T varies substantially when 9 jumps. One 
might wonder if the results here are particular to our model. In fact, brisk variations of the phase are widely shared 
by many types of models^^-'^^. 

The third case corresponds to kpa not close to m:/2. In this case, for large Xo, Q is almost zero except in the 
vicinity of a value t]c, which is very near ^/2Xo', numerically it seems that rjc is always a little less than this value. The 
maximum pumped charge is of order e but no longer close to one electron charge. Noise has a double peak structure 
around rjc- A rough qualitative picture of this can be seen in Eq. (jl9p . because, as soon as the integration contour 
does not get close to the point X{t) = Y{t) = 0, for any t, the integrands in Eqs. (|15p and ([16]) are very small. An 
example is shown in Fig 7. 




FIG. 7: Same as Fig. 5, except that kpa — 0.4. The solid horizontal line of ordinate 1 and the thick vertical line at abscissa 
Xo\/2 are just guides to the eye. 
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IV. CONCLUSION 



We have studied the influence of weak electron-electron interactions on pumped charge and noise in the adiabatic 
regime, in a mesoscopic one-dimensional disordered wire. Within the two-delta potentials barrier model, analytical 
results were obtained for the charge and noise. Results were analyzed numerically for local pumping fields with a 
harmonic dependence. Without any voltage ofi'set, at weak pumping amplitudes, interactions tend to enhance the 
pumped charge, as where / is the interaction parameter. For fairly large pumping amplitudes, it is exactly the 

reverse, Q and Qq both decrease as 77"^, but Q remains smaller than Qq by a factor P"'. At very large pumping 
amplitudes, Q and Qo are practically the same. As to the pumping noise, at weak amplitudes, it increases with 
interactions, but in the same way as the charge, so that, the Fano factor remains constant, independent of the 
interactions. For moderate pumping amplitudes, noise has a double peak structure around the maximum of pumped 
charge. For large amplitudes, the noise decreases slower than the charge, as rj^^, and for very large amplitudes, noise 
with and without interactions become approximately the same. 

As emphasised in Ref. ^8,, interactions tend to make resonances sharper, which is conducive to obtaining an almost 
quantized pumped charge. However, it is not sufficient to enclose a resonance, it is also necessary that the pumping 
contour does not go too far from the resonance. Otherwise, the noise is appreciable and the signal Q can even be very 
small. 

In the case of constant offset Xq, the behavior depends if we are close to a resonance, kpa = mT/2 in the two-delta 
potentials model. Close to a resonance, there is a region of almost quantized pumped charge where the noise seems to 
follow a Q{e — Q) behavior, reminiscent of the noise reduction in quantum wires for good transmission by a T{1 — T) 
factor, where T is the modulus of the energy transmision coefficient. Quite generally, noise breaks up into pieces due 
to fluctuations of T and those due to fluctuations of 9. We believe that in the quantized region, the fluctuations of 9 
are predominant and give rise to the Q(e — Q) behaviour. 

Interactions do help in having a charge closer to e and to reduce the noise. However, it does not change the range 
of pumping amplitudes, where quantized charge is observed, i.e. the width of the quantized region practically does 
not depend on the interactions. Qualitative arguments seem to indicate that this width scales as sin{2kpa)~^ close 
to kpd = mT/2. Further from the resonance, the maximum charge which can be pumped becomes of order, but less 
than e. Moreover, the region of quantized charge shrinks to a narrow window of pumping amplitudes around a value 
close to Xf)^/2. 

In summary, our study of noise shows that interactions tend to increase the quality of pumping. However, two 
conditions need to be met; first, to operate at certain wavevectors favouring sharp resonances and second to have a 
pumping contour which encircles the resonant point, passing not very far from it. Otherwise, only noise is produced 
and the quantization of the charge is not achieved. These restrictions were not pointed out in previous works. In 
addition, in the quantized charge region, noise vanishes as e — Q. 
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APPENDIX A 



In this appendix, we consider the limit of large pumping. In order to explain why, for offset Aq — 0, Qo behaves 
like ■q^'^ at large -q, we use Eqs. ([15]), p^ . and (|^ . Let us look first at the terms involving dX/dt. 

X, Y, as well as their time derivatives will be of order 77, except at isolated particular points. The term 

dlnsf^/dX can be expanded in powers of l/X and l/Y: 



9 In 3^2' 



dX 



cotan{2kpa) 3 cotan^{2kpa) — 1 1 + cotan'^{2kpa) cotan^(2kpa) — 1 



X 



X 



4A 



— 2 — 2 

4A Y 



2X Y 



0(77-5). (Al) 



All terms multiplied by dX /dt and then, integrated over one period give zero. Note that a term like 1/(A Y), which 
is in 77-5, would not give 0. Thus, for large 77, the term proportionnal to dX /dt in Q, behaves at least as 77^^. 

For the term involving dY/dt, the situation is simpler. dY /dt goes as 7y, but since D goes as rf' , Tq goes as f]^^, 
thus this term is at least in 77-"^. 
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Now the remainder of contributions to Q — Qq behave at least hke ry"^ for large ij; it can be seen from ((T6)) . For large 
?7, since Tq goes as ?7~^, the quantity To/(l + To(A^" — 1)) is practically equivalent to Tq ~ ?7~^. /mj^lns^^V^^} 
goes as r]~^ for large ry. The same holds for Im{dlns^i2 /9Y}. In the integral in the r.h.s. of /mj^lns^^V^^} 
is larger, by 77, than |(91ns^2V^^)P- Then, /mj^lns^^V^^} of order 77^^ whereas To is of order 77"^. Finally, the 
integrand is at least of order Im{dlnsf^ /dX}{dX/dt)To, or Im{d\iis'i!^ /dY}{dY /dt)To, which are both at least of 
order 77"^ x 77 x rj^"^ ~ ??~^- 

We now evaluate the behavior of the noise for large 77. When u ^ ep, only the low order Fourier components of 
T]{t) are important, ci and £3 will be much smaller than ep. At ci — e'2 — 0, s^'J'' ~ — e^^**''^''(l + i/X + 0(77^^)) and 
s'22 is the same but X is replaced by Y. s^^"* = — (l/2Xy) (l + i cotan{2kFa)) + 0(7/^"^). As a result, in Eq. (j23|) . 
s(ei, t) Vzs(e2, i) = Uz +0(77"^), the same holds for s(e2, t') Vzs(ei, t'), so that the trace is 0(77"^), which yields S of 
order 77"^, at least. 



APPENDIX B 



In this appendix, we show the result of calculations for the integrals Ii and I2, using Eqs. ([57)1 and ([55)1 for 
and T(x). The integrand of Ji and /2 are non-zero only if at least one of x or x' is within distance / from an Xi, 
i — 1,2. We shall need the integrals 



Mo 

K2 
Li 



exp^— ^) dx ~ 0.603, 
J exp(^ ^^^ J (1 - x)'^ da; - 1.207. 

(•1 (•! r 2 2 "I 2 

{x-y)~^dxdy ~ 7.997, 

cos|^ — (x — j/)^ I" (a; — J/)^^ da^ dy — 6.591 

Then, plugging these values into ([55)1 and ([55]) . 

Qo - e(l-TiCi) 



) 

{x~y)-^dxdy ~ 2.088, 



- 1 



-exp(7T [3 — 7^:;23T 



/I a;^ 



.2 2 



1 772 



sp(^ 7) + exp(^ -) 

— 1 7/^ — 1 . 



5o - C2(^) (e-go) [e-C3(e-Qo) 



with 



Ci = Afo ~ 0.603, 



C2 = 4^ (8X1 + 2ii)/Afo ~ 8.613, 

C3 = _J:lz^^ ^ 0.58. 



(Bl) 
(B2) 
(B3) 

(B4) 

(B5) 

(B6) 
(B7) 

(B8) 
(B9) 

(BIO) 
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APPENDIX C 



In this appendix, we show the equivalence of Eq. ([29)1 with Eqs. (24a), (24b) and (24c) of Ref. [T^. There, the noise 
power Pap between leads a and (3 was given by 



p2 ^ 

= 2-^9fi-CX(M), 

9=1 

^symfj-,-. _ Cap,q{E) + Cap,~q{E) 
^a(3,q\^) - 2 ' 

Co^UE) = 5]E[s:^(i?)s„5(i?)]Js^,(£;)s;3^(i?)]_^, 

7 S 



(CI) 

(C2) 
(C3) 



where [A\q denotes the Fourier transform at frequency quj of a time dependent quantity A. In our case, there arc only 
two leads so that indices 7 and 5 are either 1 or 2. We are interested in Pn, so we make a = /? = 1. Then, inserting 
the value of the scattering matrix elements according to Eq. leads to 
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Lg>l ^ 

On the other hand, using Eq.f^H), we obtain 

Ln>l ^ 

which is, apart from an overall factor 2, the same as Eq. (|C4[) . with q changed to n. 



(C4) 



(C5) 
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